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Abstract 

We investigate probability density functions of velocity differences 
at different distances r measured in a Couette-Taylor flow for a range 
of Reynolds numbers Re. There is good agreement with the predic- 
tions of a theoretical model based on non-extensive statistical mechan- 
ics (where the entropies are non-additive for independent subsystems). 
We extract the scale-dependent non-extensitivity parameter q(r, Re) 
from the laboratory data. 



1 Introduction 



Recently a generalization of the ordinary formalism of statistical mechanics, 
so-called non-extensive statistical mechanics, has gained a lot of interest [[l], 
[|. |3|]. While in ordinary statistical mechanics the Boltzmann-Gibbs entropy 
is extremized (subject to constraints), in the generalized formalism the more 
general Tsallis entropies S q are extremized. Given probabilities Pi of the 
microstates of the physical system under consideration, the Tsallis entropies 
are defined as 

= -!>?), (i) 

where the parameter q is the non-extensitivity parameter. The Tsallis en- 
tropies are convex, take on their extremum for the uniform distribution, and 
preserve the Legendre transform structure of thermodynamics. However, 
in contrast to the Boltzmann-Gibbs entropy, they are non-extensive (non- 
additive for independent subsystems) provided q ^ 1. 

Ordinary statistical mechanics is contained as a special case in the gener- 
alized formalism, since in the limit q — * 1 the Tsallis entropies reduce to the 
Boltzmann-Gibbs entropy, S\ = — J2iPi hip*- Recent work indicates that the 
non-extensive formalism with q ^ 1 describes many systems, including those 
exhibiting Levy- type anomalous diffusion 01], particles produced in collider 
experiments near the Hagedorn phase transition [||, |6|, |7| , and various turbu- 
lent systems |§ [| 0, [Tl| . A list of references on non-extensive statistical 



mechanics and its applications is given in [12|. Generally, it is not known 
how the parameter q depends on the internal properties of the system under 
consideration. 

In this paper we will apply the non-extensive formalism to turbulence. In 
[ |TT|1 , fully developed turbulent states were distinguished from spatio-temporal 
chaotic states that extremize the Tsallis entropies, and formulas were ob- 
tained for probability densities of longitudinal velocity differences measured 
at a distance r. Here we test the theoretical predictions by comparing them 
with experimental data obtained for turbulent Couette- Taylor flow for differ- 
ent Reynolds numbers Re and spatial scales r. For details on the experiment, 
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We will provide evidence that a slightly generalized version of the theory 
described in |TT[ yields a good fit to the experimentally measured proba- 
bility densities. These fits were achieved by varying only one quantity, the 
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non-extensitivity parameter q. We will also present the first systematic ex- 
perimental results showing how q depends on Re and r. 

The probability densities obtained from non-extensive statistical mechan- 
ics asymptotically decay with a power law with a rather large exponent w. 
In our model w is related to q by w = (4 — 2q)/(l — q). We will provide 
experimental evidence that there is a simple scaling law for the function 
w(r, Re). 

This paper is organized as follows. Section 2 describes the experiment. 
Section 3 summarizes the theory presented in fly]], and the formalism is 
slightly generalized. Section 4 compares the experimentally measured prob- 
ability densities with the theoretical predictions. Section 5 examines the 
dependence of the non-extensitivity parameter q and the exponent w on 
Reynolds number and the spatial scale. 



2 Experiment 

The experiments were conducted on a concentric cylinder system with the 
inner cylinder rotating and the outer cylinder at rest |Tj|. This Couette- 
Taylor system had a radius ratio of 0.724. Measurements were made for 
Reynolds numbers Re up to 540000, where Re = Qa(b — a)jv (Q is the inner 
cylinder rotation rate, a and b are the inner and outer cylinder radii, and v is 
the kinematic viscosity). With Re = initially, the flow exhibits a sequence 
of bifurcations with increasing Re. The last transition that has been observed 



occurs at Re = 13000 JT3[. Here we consider data for Re > 69000, where the 
flow is strongly turbulent. 

Velocity measurements were made with a hot film probe located midway 
between the two cylinders. The rms velocity fluctuations were typically only 
6% of the mean velocity; hence the data satisfy Taylor's frozen turbulence 
hypothesis, that is, temporal fluctuations recorded by the fixed velocity probe 
should accurately reflect the streamwise spatial fluctuations Jl3| . 
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Theoretical model for probability densities 
of longitudinal velocity differences 



In a perturbative approach to probability densities in fully developed 
turbulent flows was suggested, based on a small parameter a /tt, where r is 
a typical time scale of the chaotic force driving the local velocity differences 
and 7 _1 is the relaxation time to the stationary state. From the experi- 
mentally observed skewness of turbulent distributions one can estimate the 
order of magnitude of the parameter ^/yr to be about 0.1; hence a perturba- 
tive approach makes sense. Assuming that longitudinal velocity differences, 
u(r) = v{x + r) — v(x), in fully developed turbulence extremize the Tsallis 
entropies, and that large classes of chaotic relaxation processes approach the 
Gaussian limit in a universal way (see [0), one can obtain the following 



formulas 1 1 



p(u) = -L(l + /3(g-l) e ( u ))-^l (2) 
e(u) = X -u 2 - c^f L - ^u 3 ) + 0( 7 r) (3) 

where p(u) is the stationary probability density of velocity differences, and 
e(u) is a (formal) effective energy associated with the velocity difference u. 
(3 is a variance parameter that describes a (formal) inverse temperature in 
the non-extensive statistical mechanics. Z q is a normalization constant, c is 
a non-universal constant, i.e., a constant that may be different for different 
experiments and that can also depend on q. However, the functional form 
u — |m 3 of the term of order ^fyf is expected to be universal (see Jll, H). 
The parameter q depends on the distance r in an (a priori) unknown way. 

For c = and (3 = 2/(5 — 3q), the density in eq. (2) has average value 
and variance 1. However, if c ^ and j3 is still 2/(5 — 3q), then the average 
of (u) is of order 7T and the variance a = ((u 2 ) — (w) 2 ) 1 / 2 is slightly different 
from 1. An average precisely zero and a variance of unity are achieved with 
the rescaled (renormalized) distribution given by 

p(u) = ap(a(u — (u)). (5) 
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The term e(u) in eq. (|3|) stands for an effective energy in the formalism of 
non-extensive statistical mechanics. Indeed, for q — > 1, the Boltzmann factor 
p(u) ~ e~P e \ u ) is recovered. Let us here slightly generalize the approach of 
[|TT|j by considering more general effective energy levels given by 

e( u ) = ^\u\ 2a - c^rrsign{u) (\u\ a - \\u? a ^j ■ (6) 

For a = 1, eq. (|3]) is recovered, but we will allow for more general exponents 
a as well. The physical idea behind this is similar to that of the /3-model 
of turbulence |T5| , where only a certain fraction f3 of the physical volume is 



considered to contain active eddies. In the /3-model one essentially replaces 
the structure functions (|w| m ) by structure functions (\u\ am ) with a ^ 1 
related to the intermittency parameter. Similarly, in our dynamical model we 
replace \u\ by \u\ a . An exponent a slightly smaller than 1 may be interpreted 
as describing a fractal phase space and eddies that are not space-filling. The 
formalism of non-extensive statistical mechanics is designed to include such 
systems. 

In the following we will see that while the choice a = 1 yields reasonably 
good fits of the experimental data, an exponent a slightly smaller than 1 
yields the best fits to the data. 



4 Comparison with the experimental measure- 
ments 

We compare the probability density functions determined from the Couette- 
Taylor velocity measurements with the theoretical densities p(u) obtained 
from non-extensive statistical mechanics. The result is shown in Fig. 1 for 
Reynolds number Re = 540000, equivalent to a Taylor scale Reynolds num- 
ber R x = 262 |p| . The best fits were obtained by choosing a according 
to the empirical formula a = 1 — (q — 1) = 2 — q. For the strength of the 
skewness term we have chosen, c^/tt = 0.124(g — 1). Then only one inde- 
pendent parameter q is left, which is fitted for each experimentally measured 
distribution in such a way that the relative mean square deviation integrated 
over u takes a minimum. Although we only vary a single parameter q, the 
agreement with the experimentally measured densities is excellent. If a is 



4 



chosen as 1 (as originally suggested in JO]]), the agreement is still reasonable 
but not as good as for a = 2 — q. 

Other theoretical approaches to turbulent densities |Tj|-|f22|], based, for 



example, on stretched exponentials or other functional forms, usually fit only 
certain parts of the distribution (e.g., the tails), leaving other parts (e.g. the 
vicinity of the maximum) unaccounted for. Our formula yields good fits of the 
experimental data for the entire range of u- values. To demonstrate the good 
fit for all r we present both a logarithmic plot, which emphasizes the tails 
(Fig. 1(a)), and a linear plot, which emphasizes the maximum and its vicinity 
(Fig. 1(b)). To the best of our knowledge, there is no other theoretical model 
that yields fits of similar quality. Figure 2 shows that the residual difference 
between theory and experiments for the same spatial scale, r/rj « 11 (where 



rj is Kolmogorov scale ||13||), decreases with increasing Re. 

The function g(r,Re), deduced from the experiment, will be examined 
in the next section. All relevant information on the densities appears to be 
encoded in this function, which is similar to an equation of state in ordinary 
thermodynamics. 



5 Measuring g(r, Re) 

We have determined the function g(r, Re) for different spatial scales r and 
Reynolds numbers Re by minimizing the integrated relative quadratic devia- 
tion between measured distributions and theoretical distributions. For each 
Reynolds number, 38 different spatial scales were evaluated and the results 
are shown in Fig. 3. The non-extensitivity parameter q varies with both the 
spatial scale r/77 and with the Reynolds number. There is a tendency to 
smaller q for smaller Re. For small r, q approaches the value go ~ 1.185. 
For large r, q does not approach 1 but saturates at a slightly larger value, 
goo ~ 1.03, the precise value being Reynolds number dependent. Hence small 
deviations from the Gaussian distribution remain at the largest scales. This 
may be a finite-system size effect. 

The predicted probability densities for large \u\ exhibit power law decay, 
as can be deduced from eq. (H) and (|), where we neglect the term 0(y/yf) 
(our formula for this term represents a perturbative result valid only for 
\u\ < l/y/rr ~ 10): 
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p(u) ~ \u\- w (7) 

with exponent w given by 

2a 4 - 2g 



g- 1 g 



(8) 



The exponent w has the meaning that only moments (H m ) of the density 
with m < (w — 1) exist. This sounds like a severe restriction on the existence 
of structure functions, but since w is rather large (9 < w < 60 for the Couette- 
Taylor data), this effect does not contradict the experimental measurements 
of structure functions. Fig. 4 shows the function w(r, Re) determined from 
the experimental data. We observe that w exhibits simple scaling behaviour 
for medium spatial scales, 

w(r)=A^j , (9) 

where the exponent 5 depends weakly on the Reynolds number Re. We find 
5 = 0.440, 0.395, 0.360, and 0.326 for Re = 69000, 133000, 266000, and 
540000, respectively. For r/77 < 10, the experimental data saturate at w ~ 9. 
For w > 60, fluctuations become large. 

The constant in front of the experimentally determined power law is found 
to be 4.0 ± 0.1, independent of the Reynolds number. A simple argument 
for the value 4 could be as follows. Suppose we could measure the proba- 
bility distribution of ideal turbulence in an unperturbed way down to the 
Kolmogorov scale r] for infinite Reynolds number. Assuming that the scaling 
law (P) remains valid for Re — > 00 and r — > rj we obtain at r = rj the value 
w(rj) = 4, no matter what 5 is. But for turbulence to make sense at least the 
third moment (\u\ 3 ) should exist, since this is the most fundamental observ- 
able related to energy dissipation. Existence of the third moment at r — rj 
is guaranteed if w = 4 + e, where e is an arbitrarily small positive number. 
This argument suggests that the constant should be 4. 

For Re 00 the idealized small-scale turbulence at r = rj is characterized 
by the smallest possible w where it makes sense to speak about energy dissi- 
pation (the third moment). Since w = 4, the flatness factor F = (u A )/(u 2 ) 2 
as well as all higher moments at this scale would then diverge for Re 00. 



This is compatible with experimental observation [23]. 
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Figure captions 

Fig. 1. Experimentally measured probability density functions of the 
velocity differences for the Couette- Taylor experiment at Re = 540000 are 
compared with the theoretical curves p(u) (for a = 2 — q) : (a) Logarithmic 
plot, (b) Linear plot. The logarithmic plot is sensitive to the tails, while the 
linear plot is sensitive to the vicinity of the maximum. For the experimental 
curves, the distances r/rj (where r\ is the Kolmogorov length scale) are, from 
top to bottom: 11.6, 23.1, 46.2, 92.5, 208, 399, 830, and 14400. For the 
theoretical curves, the values of the non-extensivity parameter q are, from 
top to bottom: 1.168, 1.150, 1.124, 1.105, 1.084, 1.065, 1.055, and 1.038. For 
better visibility, each distribution in (a) is shifted by -1 unit along the y axis, 
and each distribution in (b) is shifted by -0.1 unit along the y axis. 

Fig. 2. Relative difference d(u) = (p t h(u) — p exp t(u)) /pth(u) between the- 
oretical and experimental probability densities at r/r] — 12. The difference 
decreases with increasing Reynolds number Re = 69000 (+), 133000 (x), 
266000 (*), and 540000 (□). For each Reynolds number the best possible fit 
is used (q = 1.148, 1.159, 1.167, and 1.168, respectively). 

Fig. 3. The fitness parameter q(r, Re) deduced from a least square fit 
of the velocity data to the theoretical probability density. The Reynolds 
numbers are 69000 (+), 133000 (x), 266000 (*), and 540000 (□). 

Fig. 4. The exponent w = (4 — 2q)/(q — 1) describing the decay rate of 
the probability density for large \u\. Scaling behavior of w(r) is observed for 
a large range of distances r. The Reynolds numbers are 69000 (+), 133000 
(x), 266000 (*), and 540000 (□). The straight lines correspond to the power 
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law w(r) = 4(r/rj) s with 5 = 0.440, 0.395, 0.360, and 0.326, respectively. 
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Fig. 1a 
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